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Abstract Thus far, there does not appear to be an agreed 
(or adequate) definition of homogeneous dark energy (DE). 
This paper seeks to define a valid, adequate homogeneity 
condition for DE. Eirstly, it is shown that as long as 7 ^ — 1, 
DE must have perturbations. It is then argued, independent 
of Wjc, that a correct definition of homogeneous DE is one 
whose density perturbation vanishes in comoving gauge: and 
hence, in the DE rest frame. Using phenomenological DE, 
the consequence of this approach is then investigated in the 
observed galaxy power spectrum - with the power spectrum 
being normalized on small scales, at the present epoch z = 0 . 
It is found that for high magnification bias, relativistic cor¬ 
rections in the galaxy power spectrum are able to distinguish 
the concordance model from both a homogeneous DE and a 
clustering DE - on super-horizon scales. 

Keywords General relativity • Dark energy • Perturbation • 
Homogeneity • Matter • Power spectrum 

1 Introduction 

Dark energy (DE) is dark, but the underlying physics of 
DE is even darker. Understanding the nature of DE remains 
a puzzle in general relativity. A long standing question is 
that: is DE actually static vacuum energy A, i.e. like in the 
standard concordance model (ACDM); or a dynamic field, 
e.g. like in the quintessence Q-lETI models (QCDM)? If 
DE is described by A then it can not have perturbations (or 
evolve) at all. Although ACDM is the best-fit model to the 
current data, other alternatives have been considered in the 
literature, e.g. a homogeneous dynamical DE. However, if 
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DE is dynamical then it can have perturbations. How then 
do we define a valid homogeneous dynamical DE? 

There does not appear to be an agreed or adequate defi¬ 
nition of homogeneous (dynamical) DE. Eor example, if the 
DE physical sound speed is = 1, then DE cannot clus¬ 
ter on sub-Hubble scales. Thus it is sometimes concluded 
that DE is approximately homogeneous (see e.g. 1^-14014 
The caveat to this assumption is that it only ensures a ‘scale- 
dependent’ homogeneity, in the sense that it makes DE ho¬ 
mogeneous only on sub-Hubble scales, but on super-Hubble 
scales DE becomes significantly inhomogeneous. This is be¬ 
cause, Csx = 1 implies that the DE density perturbations prop¬ 
agate with the speed of light; hence DE fails to cluster, and 
is perturbatively insignificant. However, on (Hubble) hori¬ 
zon scales the perturbation behaviour is different and the 
homogeneity assumption breaks down, i.e. the implicit as¬ 
sumption of no clustering in DE on super-horizon scales, 
given that = 1, is inconsistent. 

Moreover, an assumption often used for DE homogene¬ 
ity is the requirement that all its perturbations vanish m- 
Oni, i.e. by setting the DE density perturbation and velocity 
potential to (absolute) zero Q : 

4 = 0 = Vx, (1) 

where the associated evolution equations are therefore dis¬ 
carded. However, it has been pointed out that a fiuctuating, 
inhomogeneous component is the only valid way of intro¬ 
ducing an additional energy component (i.e. DE): a smooth 
(non-fiuctuating), time-varying component is unphysical - 
it violates the equivalence principle m. Moreover, it is 
known that Eq. 0 leads to a violation of the self-consistency 
of the equations of general relativity, by causing a contra¬ 
diction in the equations. Eor example, it has been shown 
that Eq. 0 leads to a false boost in the matter power spec¬ 
trum |[Il[ 6 l, and in the integrated Sachs-Wolfe effect ISlfTTIl 
(in the cosmic microwave background) - on horizon scales. 
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(Also, recently it has been shown in EB that neglecting the 
DE perturbations can lead to misleading estimation of the 
matter growth index, giving up to ^3% deviation: which is 
a significant amount, as we enter an era of precision cosmol¬ 
ogy.) Nevertheless, none of the previous works has shown 
explicitly what the inconsistency resulting from Eq. 0 is, 
nor has given any suggestions on how to solve or circum¬ 
vent this problem. 

In this work we show analytically the inconsistency re¬ 
sulting from Eq. 0 . We propose an alternative way, via the 
intrinsic entropy perturbation, to define a suitable condition 
for DE homogeneity - which corrects Eq. 0 to avoid the 
violation of the consistency of the equations of general rel¬ 
ativity, and also eliminates Csx from the equations (like in 
ACDM). It should be noted that the aim of the paper is not 
to fit the given homogeneous DE to the data (which might 
have been one other avenue to show the observational sig¬ 
nificance of the homogeneous DE), but to provide a suitable 
way to define a valid, adequate DE homogeneity condition - 
which currently seems to be non-existent. Eurthermore, for 
illustration purpose, the effects of general relativistic correc¬ 
tions (and magnification bias) in the galaxy power spectrum 
is demonstrated. 

We begin by outlining the basic background equations in 
section we give the perturbations equations in section 
In section we discuss the intrinsic entropy perturbation; 
and in section we look at DE homogeneity in general rel¬ 
ativity: discussing the ‘unphysical’ smooth DE and a ‘true’ 
homogeneous DE - illustrating their effects in the galaxy 
and matter power spectra. We conclude in section 

2 The Background Equations 

The standard acceleration equation (see e.g. 14^ ) is: 

Jf' = -ljf\l+3w), w^J^QaWa, (2) 

^ A 

where = a'/a is the comoving Hubble parameter, prime 
denotes derivative with respect to conformal time, with a 
being the scale factor; = Pa /p is the energy density pa¬ 
rameter (with over bars denoting background) of species A, 
which evolves according to the equation 

Q'^ = -3 J^{wa-w)Qa. ( 3 ) 

with pA being the background energy density of species A, 
and p is the total energy density of all the species. Similarly, 
>^A = Pa/Pa is the equation of state parameter of species A, 
which evolves by 

Wa = - 3Jif{l+WA)iclA-WA), ( 4 ) 

where c^a = PaIPa i^ adiabatic sound speed associated 
with species A, and pA is the background pressure. 


3 The General Perturbations Equations 

Here we adopt the Newtonian metric, given by 

ds^=a^ [-(l-f20)Jr]^ + (l-20)Jx^], (5) 

where 7] is the conformal time, and <P is the (Newtonian) 
gravitational potential. Note that by the choice of the met¬ 
ric ^ we assume zero (or negligible) anisotropic stress: this 
assumption is crucial for the subsequent derivations. The rel¬ 
ativistic Poisson equation is given by 

v^<P = ljr^'£nAAA, (6) 

where the comoving density perturbation Aa is given by 

Aa = 5a + ^Va = 5a- 3Je{ 1 + wa)Va, (7) 

PA 

where 5a = 5Pa/Pa and 5pA is the energy density pertur¬ 
bation. The gravitational potential is driven by the total mo¬ 
mentum density, given by 

+ = ( 8 ) 

where the 4-velocities are given by Ea 

(1 - <P, d‘VA), = a"' (1 - 4>, dW ), (9) 

with being the total 4-velocity and V is the total velocity 
potential, given by 

V = -^ynA{l+WA)VA. (10) 

We consider all species as perfect fiuids. Thus for the 
species A, the perturbed energy-momentum tensor is 

57 ’/'' = {5pA + 5pA)u'^ul + 5pAg^''+PA5g^'' 

+ {Pa+Pa)[5u’^uI + u^5uI\, ( 11 ) 

where 5pA, 5u^ and 5g^'' are the perturbations in the pres- 
sure, 4-velocity and the metric tensor, respectively. The con¬ 
servation of energy and momentum implies that 

^hY.5T^^ =0 = V^5T^\ ( 12 ) 

A 

where the second equality follows from the assumption that 
the individual fluid species do not interact directly with one 
another: they only interact (indirectly) gravitationally via the 
Poisson equation ([^. 

Thus given Eq. ( p^ the velocity potential Va and the 
comoving overdensity Aa 0 evolve by 

V^+jeVA = -<P-x^AA, (13) 

1 +Wa 

4 - 3JTwaZIa = - (1 + wa) Wa, (14) 

where we have defined the parameter Aa by 

jriA = ^^2(1 +Wa) ^ ilB{l+WB){VA-VB). (15) 

^ B^A 

The index B runs through the entire species, for a given 
(fixed) value of A. 
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4 The Intrinsic Entropy Perturbation 


The entropy of a given (thermodynamic) system or fluid, 
measures the degree of ‘disorderliness’ of that fluid; hence 
is a perturbed quantity. The intrinsic (or inherent) entropy 
perturbation of A, may be given by ll44l - ll48]| 


Pa8sa = Pa 



SPa \ 
Pa ) 


(16) 


i.e. the entropy perturbation quantifles the part of the (ef¬ 
fective) pressure perturbation that is not simply related to 
the (effective) energy density perturbation. Then the physi¬ 
cal sound speed of species A, is deflned in the rest frame 
(“rf’) of A - given by 


= 


8pA 

SpA 


> 0 , 

rf 


(17) 


where we note that this is essentially the speed of propaga¬ 
tion of the pressure perturbation 5 pa relative to the density 
perturbation 5pA- when A is at rest. Then by changing from 
some arbitrary frame into the rest frame v^|rf, given by 
x^y = x^ , this leads to a gauge transformation 
of the energy-momentum tensor: 


T^y ^ T^y\ri = T^y-^pAy^ ( 18 ) 

where the Lie derivative , with respect to the transforma¬ 
tion 4-vector is given by 

= ^"daTj^v + T^ady^"-T«yda^^ (19) 

with <C 1. Thus the (O-O)th and the (/- 7 )th components 
of the transformation yield, respectively 

spAy = spA-^^py 5pAy=5pA-^^py ( 20 ) 

and the (/-0)th or (0- 7 )th component yields 

VA|rf = V:4 + ^‘’, ( 21 ) 

with the velocity potential Va being given in the Newtonian 
gauge - where it is automatically gauge-invariant. (One ad¬ 
vantage of using the conformal Newtonian metric is that the 
resulting gravitational perturbations are automatically gauge- 
invariant.) However, when A is at rest, we have 


T^j\^ = 0=Tio\^, ( 22 ) 

where it follows that VA|rf = 0, and thus = —Va. Putting 
this in Eq. ( |20| ), we then obtain by using Eq. ( p^ that 

PpA = cIaPPa + (CsA - cIa)Pa^a, (23) 


where Aa is given by Eq. Q. Whence we obtain the intrinsic 
entropy perturbation (EH, given by 


where Aa is gauge-invariant, and consequently so is 5s a- 
Moreover, given Eqs. 0 and (|20|)-([22|), we get 

^Pa I rf — Pa^A' (25) 

This implies that the comoving density perturbation of any 
species corresponds to the density perturbation of that species 
in its rest frame. Moreover, given Eq. ( [25] ), we have Eq. ( [24| ): 
Pa5sa = (c^A ~ ^aA)^P^y ~ which therefore implies that 
the ‘intrinsic entropy’ of any species corresponds to the ‘en¬ 
tropy perturbation in the rest frame’ of the species. 


5 General Relativity and Dark Energy Homogeneity 

In this section, we analytically discuss the inconsistency in 
the equations of general relativity - resulting from Eq. 0, 
we illuminate what the inconsistency really is: (analytically) 
describing its source/origin. We then propose a suitable way 
to deflne a valid, adequate condition for DE homogeneity in 
general relativity. 

It is known that the equations of general relativity form 
a complete and consistent system. Thus an implication of 
this is that the gravitational potential evolution equation (0 
should always reduce to the Poisson equation 0. To con¬ 
firm this, it is only sufficient to show that the Poisson equa¬ 
tion at any time solves the associated gravitational potential 
evolution equation. 

Hence by taking the time derivative of Eq. ([^, and using 
Eqs. (|^, ([^, ([7])-(p3]), we get 

+wa)V^Va, (26) 

where by applying the inverse Laplacian to both sides, we 
get Eq. ^ - as required. This way, the system of equa¬ 
tions remains complete and consistent. Nevertheless, note 
that Eq. ( [^ is obtained mainly as a result of the fact that 

(27) 

A 

where given Eq. EH, it is easy to establish Eq. 

Equation ( [27] ) is essentially the statement of a ‘consis¬ 
tency condition’ for the system of equations of general rel¬ 
ativity. This condition should always hold given any correct 
set up - within general relativity. (Note that Eq. imme¬ 
diately holds true for ACDM. However, as will be shown 
subsequently, there are situations for dynamical DE where 
this condition would (i) hold true, and (ii) not hold true.) 

However, if in any adopted framework we have that this 
condition does not hold, i.e. £a -^aAa 0 , then this will re¬ 
sult in a contradiction: where we are unable to recover the 
standard gravitational potential evolution equation (0; ef¬ 
fectively we will rather have a transformation, given by 


PaSsa = (c?^-c^^)Pa4a, 


(24) 0' ^ + 


(28) 
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where 0' is given by Eq. ([ 8 ]), and 

3^^ = (29) 

^ A 

Thus Eqs. ( |^ and ( [29] ) analytically express the ‘unwanted’ 
inconsistency (or contradiction) that will arise in the phys¬ 
ical equations of general relativity, when Eq. ( j^ fails to 
hold. It should be noted that the parameter J^Aa is physi¬ 
cal, and will contribute to the comoving density perturbation 
via Eq. ( [Td] ). 

In the following subsections, we give the particular evo¬ 
lution equations and present a new definition for a ‘true’ ho¬ 
mogeneous DE. 

5.1 The particular perturbations equations 

We assume (henceforth) the late Universe - dominated by 
DE and matter (“m”), i.e. baryons and cold dark matter. Thus 
the relativistic Poisson equation is 

V20 = ^ , (30) 

where we take care to use the correct overdensities in 
order to avoid ‘unphysical artefacts’ (see e.g. 1491 ) in the 


results. 

The gravitational potential evolves by 

= -^^2 [n^v„+n,{i+w^)v^], (31) 

and the matter perturbations evolve according to 

V;„ + JfVm = -^, (32) 

+W,){Vn,-V^) = -V^Vn,. (33) 

Similarly, the DE perturbations evolve by 

V; + « = -<P- -^A^, (34) 

-{l+>v^)v2y^, (35) 


where Eqs. ([^-([^ follow from Eqs. (|6J-([T5J. 

Thus Eq. O implies that Ajc = 0 = A^, and Eq. •Ulif 
yields: —(9/2)^^^2^(l -\-Wjc)Vm = 0 , which implies that 
either (i) = — 1, or (ii) Vm = 0. But it is already taken that 

Wx 7 ^ — 1, and Vm cannot be zero (in the given gauge). Hence 
this leads to a contradiction. Usually in the literature, the A^ 
equation is merely disregarded (while keeping 7 ^ — 1 and 
Vm 7 ^ 0 ) - this is the source of inconsistency in the general 
relativistic equations, which has rightly been reported in the 
literature (see e.g. iiiiii). 



k {Mpc -^) 


Fig. 1 Plots at z = 0. Top: The DE and the matter velocity potentials, 
Vx{k) (solid line) and Vm{k) (dashed line), respectively; and their differ¬ 
ence Vm{k) — Vx(k) (dotted line). Bottom: The parameter J^Amik), for: 
clustering QCDM (solid line), and unphysical smooth QCDM (dashed 
line). The vertical line denotes the (Hubble) horizon. 

Hence given that 7 ^ — 1, and since by Eq. Q the A'^ 
equation (and hence, A^) is discarded, we have 

^ OAAA = \.M’^OA{l+W^)ym = ^^0,nAm. (36) 

A=mA 

which eventually leads to Eq. ( [^ ; thereby defying Eq. ^Tf\ . 
Obviously by Eq. ( |3^ (and preceding explanations), unless 
Wx = — 1, Eq. ([T]) leads to a contradiction - and hence a vio¬ 
lation of the self-consistency of the equations of general rel¬ 
ativity. Thus Eq. ([T]) is wrong, and the resulting ‘smooth’ DE 
is ‘unphysicaT. However, Eq. ( [3^ reveals why the ACDM 
satisfies general relativity despite the fact that all the DE per¬ 
turbations therein become zero: Wx = — 1 in ACDM. Basi¬ 
cally, provided Wx 7 ^ — 1, DE must cluster. 

However (for completeness), if we consider a (generic) 
clustering DE, i.e. with 14 7 ^ 0 and Ax 7 ^ 0 7 ^ 5x, we get 
^2mAm + = 0. It is a straightforward thing to show 

that, by taking the time derivative of Eq. and apply¬ 
ing the appropriate equations, we obtain Eq. Hence, a 
clustering DE rightly upholds the consistency of the system 
of equations of general relativity. But as previously men¬ 
tioned, J^fArn (for example) is physical and will contribute 
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to the matter density perturbation via Eq. ( [33] ). Gener¬ 
ally, for the clustering DE, the growth of on super¬ 

horizon scales will be restrained by the relative velocity po¬ 
tential, Vm — Vx, while for the (unphysical) smooth DE this 
term grows almost linearly, driven by (see Eq. ([3^). 

We illustrate these behaviours at the present epoch in 
Eig. using QCDM. These behaviours explain the matter 
power spectrum reported in ID, where it is shown that a 
smooth QCDM specified by Eq. Q, leads to a false (un¬ 
physical) amplification of the linear matter power spectrum 
on super-horizon scales. Note that, although the effect of Eq. 
0 may not appear to be significant on sub-horizon scales, 
it is nevertheless crucial to use the correct and consistent 
general relativistic equations for (generally) valid analyses. 

Henceforth, we reserve the name ‘smooth’ for the un¬ 
physical DE, defined by Eq. ([T]). 


5.2 True homogeneous dark energy 

A physical, consistent homogeneous DE should be one such 
that it maintains the consistency of the equations of general 
relativity, by upholding Eq. irrespective of the nature 
of its equation of state parameter 

Given that the entropy (perturbation) of any fluid mea¬ 
sures the degree of disorderliness in the fluid, then homo¬ 
geneity or inhomogeneity of the fluid may suitably be de¬ 
fined with respect to its entropy. Particularly, that the (net) 
intrinsic entropy perturbation of the fluid vanishes. By the 
vanishing of the intrinsic entropy perturbation, it implies 
that the net internal distortion (i.e. the total change in the 
inherent distortions - quantified by the brackets in Eq. ( [T^ ) 
of the fluid becomes zero. This way, the fluid may be thought 
to be constituted by an even distribution of equi-amplitude 
distortions; hence the fluid is homogeneous (or uniform). 

Therefore, here ‘homogeneity’ refers to ‘uniformity’, so 
that a ‘true’ homogeneous DE is not one entirely devoid 
of perturbations, but one made up of uniformly distributed 
(equi-amplitude) perturbations. Thus when the DE intrinsic 
entropy perturbation vanishes, i.e. = 0, Eq. ( [24| ) yields 

{cl-cl,)A, = 0, (37) 

where either (7) or (2) A^ = 0. It is important to 

note that the definition ( [Tt] ) is independent of the choice of 
the spacetime gauge. It may also be pointed out that, if ini¬ 
tially by a priori assumptions then automatically 

5sx = ^\ however the converse is not necessarily true: if ini¬ 
tially by a priori assumptions = 0 then it may or may 
not mean that since it can also mean that = 0 

instead. 

In what that follows, we investigate the two cases (7) and 
( 2 ), given above. 


Case 1: c 


2 

ax 


= C 


2 

5'X 


If the adiabatic sound speed is equal to the physical sound 
speed, i.e. then Eq. 0 implies that 




C.v = Wv- 


3 . 2^(1 + Wx) 


= Gx > 0. 


(38) 


where this means ?>^Wx > w^/(l + w^), and either: 

(i) Wx> —I and < 0 , or 

(ii) Wjc < — 1 and > 0 . 

Thus, unless > 0, cannot be an absolute constant. It 
may only be asymptotic to a fix value, such that ^ 0 , oth¬ 
erwise Csx becomes imaginary - and small perturbations be¬ 
come unstable. Hence, the given homogeneous DE does not 
admit = constant < 0 (i.e. negative constants). Moreover, 
conditions (i) and (ii) above, imply that can not oscillate: 
it may only be either monotonically decreasing (w^ < 0 ) or 
monotonically increasing (w^ > 0). Thus, Case 1 ( [38] ) essen¬ 
tially ‘fixes’ the DE background evolutions. 

To illustrate Case 1 ( [3^ , we consider the well known 
Chevallier-Polarski-Linder (CPE) parametrization 15011511 : 


Wx{a) = wo + Wa{l-a), (39) 

where the scale factor a = (1 +z)“\ with z being the red- 
shift; Wo and w^ are (free) constants. We consider two sce¬ 
narios of Wx'. a generic clustering DE (CPE) and a homoge¬ 
neous DE (homCPh), given by 

CPE : w;, > -1, w; > 0, A;, ^ 0, (40) 

homCPh: w^ > — 1, w^ < 0, A^ 7 ^ 0, (41) 

where we choose wq = —0.8 and w^ = —0.2 for CPE; wq = 
— 0.8 and w^ = 0.6 for homCVC, i.e. here we consider only 
the scenario (i) of Case 1 ( [38] ). We show the behaviour of w^ 
in Eig. [^ (top panel), for both the CPE ( [40] ) and the homCPL 
<ED’ respectively. 

Throughout this work, we initialize evolutions at the de¬ 
coupling epoch, given by 1 +zj = 10^. We normalize all 
the power spectra on small scales, at z = 0 : i.e. by choosing 
the same matter density parameter QrnO = 0-24 and Hubble 
constant Hq = 73 kms“^Mpc“^ for all cases. Thus all of the 
power spectra match each other at z = 0 , on small scales. 
(The advantage of this is that any clustering or GR effects 
become isolated on large scales.) We used adiabatic initial 
conditions (see Appendix [^ for the perturbations. 

We show in Eig. [^(bottom panel), the radial galaxy power 
spectrum (see Appendix^ with galaxy bias Il52l-t57l 
b = I and magnification biasEH = 1; and the matter 
power spectrum at the present epoch, i.e. z = 0. By our 
normalization, we see that can be approximated on sub- 
Horizon scales by q^Pm, with q = ^/2A. Moreover, although 
DE clusters in the CPE, we see that it rather leads to higher 
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k (Mpc ^ ) 

Fig. 2 Top panel: The evolution of the equation of state parameter 
for the CPL ( [^ (solid black) and the homCVT ( |4T) (dashed black). 
Bottom panel (at z = 0): The radial (i.e. p = 1) galaxy power spectrum 
(solid lines) with galaxy bias b= I, magnification bias = 1; and 
the matter power spectrum - given by q^Pm (dashed lines); q = 


power on horizon scales in both and Pm, i.e. relative to 
the homCPh for equal values of (see II121I131I221I481I58II - 
(781, for the effects of c^^). This may be owing to the be¬ 
haviour of Wx for homCPh, which suggests that DE sets in 
relatively earlier for the homCPh - hence causing the matter 
perturbations to have less time to cluster, thereby resulting 
in relatively lower power spectra. Obviously, Case 1 
implies that for equal values of the difference between 
a homogeneous DE and a clustering DE is mainly governed 
by the background, with little to do with the perturbations. 
However, one may expect that this difference strongly per¬ 
tains perturbations, and that a homogeneous DE results in 
higher power spectra on large scales, relative to a clustering 
DE - given that the perturbative effect of the homogeneous 
DE should be negligible (or even absent). 

We also observe the dependence of homCPh on i.e. in 
the power spectra, with smaller values of resulting in 
more power suppression - since in which case the DE per¬ 
turbations are able to cluster earlier and on smaller scales; 


thus suppressing most of the matter growth. However, we 
see the effect of the general relativistic (GR) corrections (H 
|52l-|l52l, CD -(931 in they lead to a sizeable power 
boost (relative to Pm) on horizon scales. Moreover, we ob¬ 
serve that the GR corrections also result in significant dif¬ 
ferentiation of the given DE scenarios. 

Nevertheless, for self-consistent models, e.g. the QCDM 
(specified by a scalar field (p) which evolve along a potential 
given (generically) by U{(p) ^ constant, we have 


1 + 


/ 2 ^ 
3 ^ 9 ' ^ 


(42) 


where f/|(p = dU(o)/d(p, with and being as defined 
in sections [^and|^ respectively. Thus by Eq. ( |42| ), ^ 

which then disallows Case 1 ( [38] ). Therefore, may os¬ 
cillate (see, e.g. (41I161I94II ) or take any behaviour. However, 
one may choose to fix f/( 9 ) = constant, thereby satisfying 
Case 1 ( [38] ) (in principle). Practically though, this choice 
leads to violating Case 1 ( [^ , i.e. by becoming Wx = — I 
for 0 < z < 100 , in which case the perturbations equations 
become unsolvable numerically. Thus Case 1 ( [3^ is ‘im¬ 
practical’ for the QCDM. 

In general. Case 1 ( [38] ) is unsatisfactory, given that it still 
depends on the behaviour or choice of c^^>0. 


Case 2: Ax = 0 

On the other hand, if the DE comoving density perturbation 
vanishes, i.e. Ax = 5x — 3^(1 -i-Wx)Vx = 0 , then it implies 
that the correction to Eq. ([T]) is suitably given by 

5x = 3Jf{l^Wx)Vx, Vx^O, (43) 

where consequently, A^ = 0. Thus Eq. ( [4^ implies that the 
homogeneous DE should not have any density perturbations 
Ax in comoving gauge, but may posses the fractional den¬ 
sity fiuctuations 5x which generate peculiar velocities with 
potentials Vx. 

Note that given Eq. Case 2 ( [43] ) implies that a ho¬ 
mogeneous DE has zero density perturbation in its rest frame, 
i.e. 5px\rf = pxAx = 0. This is the physical statement of Case 2 
which is a crucial statement - as it suggests that DE 
homogeneity or inhomogeneity should be defined relative to 
the DE rest frame, i.e. by whether or not its density pertur¬ 
bations vanish in its rest frame. 

In fact. Case 2 ( [4^ readily holds for the QCDM when 
5sx = 0 (i.e. from Eq. ([T^); in which case we have 

^ax^PfP — ^P(p — ^P(p -\-Wx){c^x~ (44) 

where by collecting terms with 5p(p to one side and divid- 
ing through by - l)p,p, we get 4 = 3^(1 + We 



























Dark energy homogeneity in general relativity: Are we applying it correctly? 


7 


have used Eq. ( |4^ ; 9 '^ = a^{l -^Wx)p(p, = P(plp(p with 

the perturbations Vx = —5(p/(p\ 5x = Sp(p/p(p and 

Sp(p = 5p(p — 2U^(p5(p, (45) 

dp^ = {d(p' - 9'4>') + Ui^5(p. (46) 

Besides, by considering the definition of the intrinsic en¬ 
tropy perturbation specifically for QCDM given by F iflTt . 
which relates to the entropy perturbation of section by 
r = Wx 5sx/(c^x ~ ^ax) = - where the second equality fol¬ 
lows from Eq. ( [^ , then it automatically follows that Ax = 0 
when we set F = 0 (which is given therein as the adiabaticity 
condition for quintessence). Thus a clustering quintessence 
will have Axy^O and a homogeneous quintessence will have 
Ax = 0: both scenarios having the same (background) equa¬ 
tion of state parameter — 2a^U)l{(p'^ + 2a^U), 

which varies (in general) by — 1 <Wx< 1 . 

It should be pointed out that the discussion on Eqs. (|44|)- 
( |46| ) is not implying that QCDM models are generically in 
the form of Case 2 ( [43] ), nor are we claiming that there is any 
other particular DE model (or class of models) that naturally 
exists in such form. Instead, Case 2 ( |4^ is rather a propo¬ 
sition for a ‘general’ homogeneity condition which may be 
applied to any (dynamical) DE model. (Notice that Case 2 
( [43] ) was arrived at - mainly via Eq. ( [Tt] ) - without assuming 
any DE models). Hence the QCDM models are only used as 
a reference or an example of a particular, well known DE 
model that Case 2 ( [43] ) may easily be applied to. Moreover, 
Case 2 ( [43] ) applies to the CPE parametrization - which ap¬ 
proximates (canonical) scalar-field DE models. 

In general. Case 2 ( |4^ suggests that DE may be ho¬ 
mogeneous only when the observer is comoving with the 
source, irrespective of the DE background specifications. 
Thus a homogeneous DE may have perturbations, but only 
such that these perturbations combine to cancel out in co¬ 
moving gauge (and hence, in the DE rest frame). Moreover, 
the matter power spectrum physically makes sense only when 
computed in comoving gauge, since otherwise, it becomes 
gauge-dependent and varies with the observers on large scales 
(see e.g. ima). Thus Case 2 being defined in co¬ 
moving gauge, can lead to (physical) observable implica¬ 
tions in the power spectrum. Moreover, the effect of Case 2 
( |4^ on the matter perturbations will be imposed directly, 
rather than indirectly via the background evolutions - as in 
Case 1 ( [^ . This way, the imprint of the given homogeneous 
DE will bear directly on the growth of structure. 

An important advantage of Case 2 over Case 1 ( [38| ) 
is that, unlike Case 1 ( [3^ , Case 2 ( [43] ) permits an arbitrary 
background behaviour for the given DE: may take any 

nature (constant or otherwise). This is important for mod¬ 
els with either oscillatory or constant A further advan¬ 
tage of Case 2 ( [43] ) is that it eliminates the dependence of 
the perturbations on (via Eq. ([34|)), given that = 0 or 
5x = { \ ^ Wx)yx- Hence the given homogeneous DE is 



Fig. 3 Plots at z = 0: The radial (i.e. jd = 1) galaxy power spectrum 
Pg (solid lines) with galaxy bias b=l and magnification bias = 1; 
and the matter power spectrum - given by q^Pm (dashed lines) with 
q = y/Fl. The vertical line denotes the horizon at the given redshift. 

completely independent of the choice or nature of (just 
like in ACDM). This further removes the risk that accom¬ 
panies a bad choice or wrong modelling of Moreover, 
it also reduces the parameter space that needs to be con¬ 
strained. 

To illustrate Case 2 ( [43] ), we use only the generalized 
phenomenological model ( [39] ) (ignoring the particular sce¬ 
nario of the QCDM). Hereafter we denote clustering DE by 
ACDM - given by Eq. ( [40] ), and denote the associated ho¬ 
mogeneous DE by homXCDM - which has the same back¬ 
ground parameters as ACDM, but with its perturbations be¬ 
ing specified by Case 2 ( |4^ . We use Csx = 1 for all numerical 
computations (i.e. for ACDM). 

We show in Eig.j^the galaxy power spectrum with 
b = I and ^ = I, and the associated matter power spec¬ 
trum Pryi '. for ACDM and homXCDM, at z = 0. We see that 
on sub-horizon scales, P^^^ can be approximated by q^Pm, 
where q = v/zT. Moreover, unlike the results by Case 1 ( [38] ) 
(see Eig. [^, where the clustering DE results in large-scale 
boost in the power spectra relative to the homogeneous DE, 
here we see that although = 1 for the ACDM, we get 
large-scale power suppression in both P^^^ and P^ relative 
to those for homXCDM (and ACDM) - i.e. the power spec¬ 
tra for ACDM are lower than those for homXCDM, on hori¬ 
zon scales. This implies that on horizon scales, the effect of 

in ACDM is less significant and hence the DE perturba¬ 
tions are able to cluster enough to suppress the growth of the 
matter perturbations. On the other hand, the DE density per¬ 
turbations vanish on all scales for homXCDM (in comoving 
gauge); thus the matter perturbations are able to grow more. 
Consequently, we get the relative boost in the power spectra 
in homXCDM. 

Moreover, we see that the ACDM gives a sizeable devia¬ 
tion in relative to homXCDM, on super-horizon scales. 
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(Note however that in reality, this statement depends on: (1) 
the cosmic variance on these scales, and (2) the error bars 
achievable by a given survey experiment. But for the pur¬ 
pose of this work, we leave out (throughout this work) any 
exact experimental aspects.) This deviation illuminates the 
sensitivity of GR corrections to changes in this sensi¬ 
tivity will be crucial in discriminating the ACDM from a 
dynamical homogeneous DE model - with the future large 
scale surveys. We also observe that. Pm of homXCDM is 
identical to that of ACDM on all scales. This reveals that, 
the linear matter power spectrum is incapable of distinguish¬ 
ing a dynamical homogeneous DE (given by Case 2 
from ACDM, on large scales - when their power spectra are 
normalized on small scales (at the given epoch). 

In the intensity mapping of neutral hydrogen (HI), usu¬ 
ally the individual radio sources are not counted: only the 
diffuse 21 cm line emission of a number of sources is de¬ 
tected (sol. Then given that the volume distortion in the 
observed density perturbation mainly leads to the ‘amplifi¬ 
cation’ of the number of sources, it is therefore taken that 
the volume distortion does not contribute in HI intensity 
mapping: hence, the magnification bias is often set to = 
1, which results in the volume distortion being cancelled 
out in the observed density perturbation. (Pure galaxy num¬ 
ber count surveys - in which mainly individual galaxies are 
counted - correspond to ^ = 0: thus eliminating the cosmic 
magnification of the galaxies.) 

In Eig.j^we illustrate the effect of the magnification bias 

We show the ratios of the radial (i.e. with /i = 1) galaxy 
power spectrum at the epoch z = 0.1, with galaxy bias 
/? = 1: for — 1, 0.5, 0.9, 1. In the top panel, we give 
the ratios for ACDM relative to ACDM. Obviously, we see 
that the ratio of the two models varies with different values 
of the magnification bias - which is not surprising since the 
relation between the magnification bias and the power spec¬ 
trum is non-trivial (see Eqs.([52|)-([54|)). If, for example, ^ 
had a simple scaling relation to i.e. if ^ appeared in 
pobs as a multiplicative factor, then the ratio between 

the given models will remain the same for all values of 
Moreover, the matter density perturbation and velocity po¬ 
tential, and the gravitational potential (hence, the GR correc¬ 
tions) in ACDM will inherently differ from those in ACDM: 
owing to the difference in the DE perturbations, and the DE 
equation of state parameter. At low z, the matter overdensity 
will be larger in ACDM than in ACDM on horizon scales. 
Consequently, for a given the amplitude of P^^ is larger 
in ACDM than in ACDM - thus the ratio of the two models 
is less than unity (on horizon scales), as seen in the plots. 

We see that the given models start to differentiate from 
each other on horizon scales, whereas on scales well within 
the horizon they match each other - which is mainly owing 
to our normalization (at z = 0). Eurthermore, we observe that 
P^^ in ACDM is consistently suppressed relative to that of 


< 








k (Mpc ^) 


Fig. 4 The plots of the ratios of (with jU = 1) at z = 0.1, with 
galaxy bias b = \. The panels show the ratios for: ZCDM relative 
to ACDM (top), homXCT>M relative to ACDM (middle) and XCDM 
relative to homXCDM (bottom). The lines denote: ^ = —\ (solid), 
^ = 0.5 (dashed), ^ = 0.9 (dot-dashed) and ^=l (dotted). The ver¬ 
tical line denotes the horizon at the given redshift. 


ACDM as the magnification bias increases: from ^ = — 1 to 
^ = 1. This implies that the value of the magnification bias 
can be crucial in differentiating a clustering DE from the 
cosmological constant (at the given z), given that the ratio 
between the models is sensitive to the values of More¬ 
over, note that for a given value of the large-scale ampli¬ 
tude of the power spectrum will be higher in ACDM than in 
ACDM (as previously explained): A has no perturbations to 
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suppress the matter perturbations; moreover, the equation of 
state parameter of A is weaker than (for dynamical DE) 
- noting that the bigger the equation of state parameter, the 
stronger the DE. 

Eig. 1^ (middle panel) also shows the ratio of for 
homXCDM relative to that of ACDM: for the given values 
of The results are similar to those between ACDM and 
ACDM (in the top panel), except that the amplitude of the 
ratios in the middle panel are relatively higher - for a given 
value of This is mainly a clustering effect of DE: unlike 
in homXCDM, the large-scale clustering of DE in ACDM 
leads to the suppression of the matter perturbations (and 
thus, GR corrections) - for the same hence resulting in 
a relatively lower galaxy power prediction in the ACDM, on 
the given scales. Moreover, just as in the top panel, the ratios 
of for homXCDM relative to ACDM also suggest that 
the value of ^ can be crucial in differentiating a homoge¬ 
neous DE from the cosmological constant. 

The bottom panel of Eig.H shows the ratios of for 
ACDM relative to horuXCDM. Similarly, we see that as ^ 
increases, ACDM becomes consistently suppressed relative 
to homXCT>M, i.e. the bigger the value of the magnification 
bias, the wider the difference between the clustering DE and 
the homogeneous DE (and vice versa). This suggests that fu¬ 
ture surveys that depend on cosmic magnification - e.g. the 
HI intensity mapping (see e.g. nil 8 ()l[ 88 l . which correspond 
to = 1 ) - will be useful in distinguishing or identifying 
a homogeneous DE from a clustering DE in the large scale 
analysis, particularly on horizon scales. Moreover, we ob¬ 
serve that the ratios (in the three panels) grow with decreas¬ 
ing suggesting that at a particular z, GR effects become 
enhanced as cosmic magnification bias decreases. 

Note that, as previously mentioned, the various ratios in 
the top and the middle panels (Eig. |^, respectively, inher¬ 
ently contain the effect of the respective equation of state 
parameters of the given models - i.e. for a value of the 
difference between successive ratios is not only owing to the 
effect of the perturbations, but also owing to the background 
difference of the models (via Wx). However, in the bottom 
panel, the ratios are resulting mainly from the difference in 
the perturbations of the models - in the comoving gauge. 
The imprint of is effectively factored out, given that we 
use the same for both ACDM and homXCDM\ thus re¬ 
vealing mainly the consequence of the Case 2^. 

In Eig. we repeat the plots of Eig. but here at z = 1. 
We observe that, generally (at the given z), on scales near 
the horizon the ratios maintain a consistent decrease with 
increasing However, in the top and the middle panels, 
respectively, the ratios grow on super-horizon scales: with 
the ratio for — 1 being the lowest, and that for = 1 
being the highest. This may be attributed to the magnifica¬ 
tion effect (i.e. terms proportional to in P^^)\ which in¬ 
creases on large scales at z ^ 1 , being boosted by stronger 



Fig. 5 The plots of the ratios of (with jU = l)atz=l, with galaxy 
bias b=\. Line notations are as in Fig.[^ 


Wx> — (The clustering effect of DE generally diminishes 
with increasing z - hence at high z the DE effect is mainly 
governed by the (background) equation of state parameter.) 
Thus P^^ becomes enhanced for increasing However, in 
the bottom panel the effect of is factored out, and only the 
effect of ^ is seen. Moreover, Eigs.j^ and [^generally show 
that for a given magnification bias, the ratios are higher on 
super-horizon scales at z = 1 than at z = 0 .1; hence implying 
higher GR effects at z ^ 1. 

In general, unlike Case 1 ( [38] ) which further depends on 
the perturbation modes being well within the horizon and 
1, Case 2 ( [43] ) is a definitive condition for DE homo¬ 
geneity: once it is chosen, no other requirements are needed. 
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Moreover, apart from solving the problem posed by Eq. O 
and the elimination of from the equations, Case 2^ 
can also conveniently admit particularly w^ — — I while still 
allowing the DE perturbations. Thus Case 2 is more ro¬ 
bust, and is hereby considered as the right DE homogeneity 
condition. 


6 Conclusion 

We have shown analytically that the DE homogeneity as¬ 
sumption = 0 = Vjc (with the evolution equations being 
discarded) violates the self-consistency of the equations of 
general relativity. We showed in Newtonian gauge that, un¬ 
less the equation of state parameter of the given homoge¬ 
neous DE is strictly = — 1, this assumption introduces a 
contradiction in the equations of general relativity. In essence, 
provided — I, DE must have perturbations. 

We have proposed a correct homogeneity condition for 
DE, which is valid irrespective of the nature of the DE equa¬ 
tion of state parameter or spacetime gauge - by supposing 
that the DE intrinsic entropy perturbation vanishes: which 
leads to the vanishing of the DE overdensity in comoving 
gauge (and hence, in the DE rest frame). Thus we correct 
the wrong homogeneity assumption (5^ = 0 = 14 ) by the fol¬ 
lowing: 5x = 3^(1 ^Wx)yx, with 14 0. A homogeneous 

DE hence, is given not as one devoid of perturbations, but 
rather as one with vanishing density perturbations in comov¬ 
ing gauge - i.e. one with zero density perturbations in its rest 
frame. Such kind of DE is not impractical. 

Using a phenomenological DE model, we investigated 
the consequence of our approach in the observed galaxy 
power spectrum. By normalizing the models at the present 
epoch on small scales, we found that: a clustering DE, a 
homogeneous (dynamical) DE and the cosmological con¬ 
stant, are each distinguishable from the others in the ob¬ 
served galaxy power spectrum, on horizon scales - suitably 
for high cosmic magnification bias. 

Moreover, the results show that for a given magnifica¬ 
tion bias, GR effects in the galaxy power spectrum become 
enhanced with increasing redshift. 


A Adiabatic initial conditions 


All evolutions in this work are initialized at the photon decoupling 
epoch, z = Zd- We use adiabatic initial conditions, which follow from 
the vanishing of the relative entropy perturbation 111171441 - 1471 . given 
atZd by 

5, = (l + w,)5„. (47) 

By using that at Zd, we have 


Vx = V^, 

then we obtain, i.e. given (( 43 , that 

= (1 + Wx)^m- 


(48) 


(49) 


These equations together with the Einstein de Sitter initial condition 
^'{zd) = 0, lead to the initial perturbations 


Vm{k) = 


-2C 

3(1+ Q.xWx) 
-2 




3(1+ flx^x) ^ 


^dik). 


(50) 

(51) 


where we take <Pd = ^(zd) as given by [T]. 


B The galaxy power spectrum 


In order to adequately account for the correct galaxy distribution on 
large scales, we use the observed galaxy density perturbation I1I52I 
I54I55I79I821 to compute that galaxy power spectrum ^^Vwhich is 
approximated in the flat-sky limit (in Fourier space) by I43l^ 



{b + fn^f + 2ib + fn^) 


^ 2 


(52) 


where Pm is the matter power spectrum 111431 : fi is the cosine of the 
angle between the line-of-sight and the wavevector k, with ^ = |k| be¬ 
ing the wavenumber; v = k/J^f is a dimensionless parameter, and 


^ = x^ 


4^-b -1 + —+ 2 ^_+. 

4^ be i+^2+^ 

h —IQ _ 9 -^_ L 


0 ' 




0 

Am 


(53) 

(54) 


where / = CVm!Am) (see also (33). Note that in ACDM, this 
/ reduces to the standard linear growth rate of matter energy density 
perturbation. We have neglected all integral terms, and assumed con¬ 
stant comoving galaxy number density - thus the galaxy evolution bias 
be = 0-,r is the comoving distance at the observed z. 
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